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ABSTRACT: The Hep process is a weak-interaction reaction, 3 He + p — + 4 He + e + + v e , 
which occurs in the sun. There is renewed interest in Hep owing to current experimental efforts 
to extract from the observed solar neutrino spectrum information on non-standard physics in 
the neutrino sector. Hep produces highest-energy solar neutrinos, although their flux is quite 
modest. This implies that the Hep neutrios can at some level influence the solar neutrino 
spectrum near its upper end. Therefore, a precise interpretation of the observed solar neutrino 
spectrum requires an accurate estimate of the Hep rate. This is an interesting but challenging 
task. We describe the difficulties involved and how the recent theoretical developments in 
nuclear physics have enabled us to largely overcome these difficulties. A historical survey of 
Hep calculations is followed by an overview of the latest developments. We compare the results 
obtained in the conventional nuclear physics approach and those obtained in a newly developed 
effective field theory approach. We also discuss the current status of the experiments relevant 
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to Hep. 

CONTENTS 

INTRODUCTION 3 

EARLIER CALCULATIONS 8 

RECENT CALCULATIONS - THEORETICAL FRAMEWORK 11 

Standard Nuclear Physics Approach (SNPA) 11 

Effective Field Theory (EFT) 13 

Nuclear EFT in the Weinberg Scheme - A-Counting 16 

Nuclear EFT in the KSW Scheme - Q-Counting 17 

Hybrid EFT 17 

EFT* or MEEFT 18 

NUMERICAL RESULTS 19 

Hep Calculation Based on SNPA 19 

Hep Calculation Based on EFT* 19 

Off-Shell Problem 20 

COMPARISON WITH EXPERIMENTAL DATA 22 

RELATED TOPICS 24 

The pp Fusion Process 24 

Neutrino- deuteron reactions 25 

Hen process 26 

SUMMARY AND OUTLOOKS 26 



2 



The Solar Hep Process 3 
1 INTRODUCTION 

The Hep process, 3 He + p — > 4 He + e + + f e , is one of the thermonuclear reactions 
that occurs in the sun. To explain why this specific process is of current interest, 
we first briefly describe the standard solar model, the solar neutrinos and neutrino 
oscillations. 

The sun generates its energy by converting four protons into an alpha particle, 
4p — > 4 He + 2e + + 2i/ e , via chains of thermonuclear reactions caused by weak, 
electromagnetic, or strong interactions. Thepp-chain, shown in Fig. 1, represents 
by far the most important scheme by which the Ap — > 4 He burning takes place in 
the sun. To establish how these reactions actually proceed in the sun, one must 
carry out a detailed simulation in which the radial profiles of the mass density, 
temperature, chemical composition, etc., are determined in such a manner that 
hydrostatic equilibrium is satisfied and the empirically known solar properties 
come out correctly. The principal inputs that go into this simulation are the 
nuclear reaction rates, equation of state, elemental abundances, and radiative 
opacity. Over the past four decades a great deal of effort has been invested 
in this subject, and out of this endeavor has emerged a quantitative model of 
the sun, called the standard solar model (SSM) (1-5). Among many quantities 
determined by SSM are the time rates of the thermo-nuclear reactions occurring 
in the sun (3); Figure 1 indicates the predicted branching ratios of the various 
paths involved in the pp-chain. Among the reactions featured in Fig. 1, five are 
weak-interaction processes that emit solar neutrinos, and SSM predicts the flux 
<p u from each source of the solar neutrinos (3); this prediction is shown in Fig. 2. 
Studying the solar neutrinos is very important for two reasons. First, it gives 
direct information about the physics of the solar interior, since the neutrinos 
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exiting the sun experience hardly any interactions with the solar medium other 

than refractive effects (related to the MSW effect to be discussed later). This 
should be contrasted with the behavior of the photons, which interact with the 
solar medium so many times that, by the time they reach the surface (after 
~40,000 years !), they do not carry much information about the solar interior. 
Second, the solar neutrinos can provide valuable information on the properties 
of the neutrinos themselves; the sun is an extremely strong neutrino source and 
hence can be highly useful for neutrino physics. 

The first measurement of the solar neutrinos was done by Davis and his col- 
laborators (6), who used a 37 C1 target. The results indicated that the sun indeed 
emits neutrinos whose flux is in approximate agreement with the SSM predic- 
tion (7), which supports the basic idea of the thermonuclear origin of the solar 
energy. At a more quantitative level, however, the measured flux was signifi- 
cantly lower than predicted by SSM. This deficit, or "solar neutrino problem", 
was also confirmed by water Cerenkov counter experiments at the Kamiokande 
(8) and Super-Kamiokande (9), by gallium-target experiments (10, 11), and by 
heavy-water Cerenkov counter experiments at the Sudbury Neutrino Observatory 
(SNO) (12). It is to be noted that because of different detection threshold ener- 
gies (see Fig. 2), these experiments are sensitive to different regions of the solar 
neutrino spectrum. If we denote by R the ratio of the observed event rate for 
a given solar neutrino detection experiment to the event rate expected from the 
SSM prediction, the current status of the solar neutrino problem is summarized 
as follows: R = 0.34 ± 0.03 for the chlorine experiment (13); R = 0.465 ± 0.015 
for the water Cerenkov counter experiment (14); R = 0.54 ± 0.03 for the gal- 
lium experiments (15-17); R = 0.35 ± 0.02 (18, 19) and R = 0.32 ± 0.02 (20) 
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for the heavy-water Cerenkov counter experiments. It should be mentioned that 

the errors attached to the above values of R only include experimental errors. 

Obviously, the degree of seriousness of the solar neutrino problem (R < 1) hinges 

on the accuracy of the SSM predictions. The latest discussion of the errors to 

be assigned to the SSM predictions (3) finds it extremely unlikely that the solar 

neutrino deficit can be attributed to the uncertainties in SSM. This conclusion is 

further corroborated by highly stringent constraints imposed by the helioseismo- 

logical data (21). 

The above discussion presupposes that the neutrinos created in the sun travel to 
the terrestrial detectors without changing their identity (or flavor). Let us recall 
that there are three distinct neutrinos, - electron neutrinos (z/ e ), muon neutrinos 
(f M ), and tau neutrinos (u T ) — and that it is the electron neutrinos that are 
produced in the sun. If there exists a mechanism (22, 23) that changes electron 
neutrinos muon neutrinos or tau neutrinos before they reach a terrestrial detector 
that detects only electron neutrinos, then there would be an effective deficit of 
solar neutrinos. This transmutation of the neutrino flavor, called neutrino os- 
cillations, signals physics that goes beyond the well-established standard model 
of particle physics, and therefore its experimental verification is of paramount 
importance. Neutrino oscillations can occur either during the neutrino's prop- 
agation in vacuum (22,23), or as the neutrinos travel in matter and experience 
refractive interactions with the medium (MSW effect) (24). 

Now, neutrinos can be detected either via charged-current (CC) reactions or via 
neutral-current (NC) reactions. Since a CC reaction involves the change v x — > x 
(where x = e~, (J,~, or t~), it can occur only for the electron neutrino; the muon 
and tau-lepton are too heavy to be created by solar neutrinos. Meanwhile, an NC 
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reaction that involves v x — ► u x occurs with the same amplitude for any neutrino 

flavor x. At SNO, the CC reaction v e d — > e~pp was used to register the electron 

neutrino flux (f> Ve , whereas the NC reaction u x d — > v x np was used to determine 

the total neutrino flux, <p u ^ = (frv e +4>v )1 +4>v T - The NC reaction data (18) showed 

that (j> v> T agrees very well with the SSM prediction (3), whereas the CC reaction 

data (12) indicated R = 0.347 ± 0.029. These results have firmly established 

flavor transmutations in the solar neutrinos. (For the evidence obtained from 

comparison of the SNO CC reaction data and the Super-Kamiokande data (14), 

see Reference (12).) Independent evidence for neutrino oscillations is known from 

the study of the atmospheric neutrinos at Super-Kamiokande (25), and from the 

study of the reactor neutrinos at the KamLAND (26). 

Neutrino oscillations occur if a neutrino state produced in a weak-interaction 
process ("weak eigenstate") is different from an eigenstate of the Hamiltonian 
(mass eigenstate), and if the mass eigenstates of different neutrinos are not de- 
generate. It is conventional to parametrize the former aspect in terms of mixing 
parameters, and the latter in terms of differences between the neutrino masses 
squared (5m 2 's). Now that the existence of neutrino oscillations has been estab- 
lished, the next important challenge is to determine the accurate values of the 
mixing parameters and 5m 2, s, quantities that should carry valuable information 
on new physics beyond the standard model. (For a recent survey of this topic, 
see e.g., Reference (27).) The great importance of this determination makes it 
highly desirable to assemble an over-constraining body of data, and this is one 
reason why detailed studies of the solar Hep process can be important. 

In discussing Hep, it is convenient to use Bahcall et aVs latest treatise on 
the SSM (3) as a basic reference (to be called BP00). According to BP00, the 
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neutrino flux due to Hep is 4> v (hep) = 9.4 x 10 3 cm^s -1 , which is seven orders of 

magnitude smaller than the pp-fusion neutrino flux, and three orders of magnitude 
smaller than the 8 B neutrino flux 4> U ( 8 B); the smallness of <j) v {hep) can also be 
seen from Fig. 2. (For the radial distribution of sites of Hep neutrino generation 
inside the sun, see Fig. 6.1 in Reference (2).) Because of its extremely small 
branching ratio (see Fig. 1), Hep in fact does not affect solar model calculations. 
So why does it interest us ? One reason is that Hep is a potential source of 
useful information on non-standard physics in the neutrino sector. Hep generates 
neutrinos having maximum energy E™ ax (hep) = 18.8 MeV, which is higher than 
that of the 8 B neutrinos, i?™ a:E ( 8 B) = 17 MeV; thus the Hep neutrinos near 
the upper end of their spectrum represent the highest-energy solar neutrinos 
(see Fig. 2). Solar neutrino detectors such as Super-Kamiokande and SNO can 
determine the spectrum of the solar neutrinos in a region dominated by the 8 B 
neutrinos. Meanwhile, the shape of </v( 8 B) is independent of solar models to an 
accuracy of 1 part in 10 5 (28). Therefore, in the absence of Hep neutrinos, any 
deviation of the observed <j) v from 4> U { 8 B) in the higher E v region reflects the non- 
standard behavior of neutrinos. Solar neutrino experiments that approach the 
level of precision needed for studying this deviation have already been reported 
from Super-Kamiokande (14,29,30) (see below). Apart from the ramifications 
for neutrino physics, the study of Hep neutrinos is also important as a possible 
additional check of the SSM itself. 

In interpreting these and future experiments, we need to know accurately to 
what extent the Hep neutrinos can affect (j) v in the 8 B neutrino region (1,31-33), 
and for this we must make a reliable estimation of the Hep cross section. This 
task, however, turns out to be extremely challenging. For one thing, although the 
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primary Hep amplitude is formally of the Gamow- Teller (GT) type, the usually 

dominant one-body GT amplitude is strongly suppressed for Hep (see below). 

Furthermore, the two-body corrections to the "leading" one-body GT term have 

opposite sign, causing a large cancellation. It is therefore necessary to calculate 

these "corrections" with great accuracy, which is a highly non-trivial task. Thus, 

from a nuclear-physics point of view, Hep presents a difficult yet very intriguing 

challenge. 

In what follows, we first present a history of Hep calculations, explaining in 
more detail the nature of the difficulties involved in Hep calculations. We then de- 
scribe how the recent theoretical developments in nuclear physics have enabled us 
to largely overcome these difficulties. After reporting the latest results obtained 
in the so-called standard nuclear physics approach, we highlight the results ob- 
tained in a newly developed effective field theory approach. We then describe the 
current status of experimental information on the Hep neutrinos. At the end we 
discuss several electroweak processes closely related to the Hep calculations. 

2 EARLIER CALCULATIONS 

An illuminating survey of the earlier Hep calculations can be found in References 
(33,34). The Hep reaction rate can be conveniently expressed in terms of the 
astrophysical S-factor defined by S(E) = Ea(E) exp(47ra/t> rc i), where a(E) is 
the Hep cross section at center-of-mass energy E, v re \ is the relative velocity 
between p and 3 He, and a is the fine structure constant; S(E) is a smooth 
function of E that remains non-vanishing as E — > 0. The first Hep calculation in 
1952 by Salpeter (35) was based on the extreme single-particle picture and only 
considered the overlap between an s-wave proton scattering wave function and 
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a Is neutron state in 4 He. This simplified treatment led to a large value for S, 

5(0) = 630 x 10~ 20 keV-b, and this value was used by Kuzmin (31), who was 

the first to discuss Hep in connection with solar neutrinos. Werntz and Brennan 

(36) pointed out the drastic suppression of the Hep rate due to a specific feature 

of the initial and final nuclear wave functions. The dominant component of 4 He 

has the orbital configuration (Is) 4 , which is totally symmetric, i.e., a state with 

[4] orbital permutation symmetry. In general, the Pauli principle dictates that a 

spin-isospin wave function accompanying an orbital function with [4] symmetry 

must be totally antisymmetric ([1111] state) and hence must have S = T = 0. 

The contraposition of this property implies that the p- 3 He state with isospin 

T = 1 cannot have [4] orbital symmetry. Meanwhile, the one-body GT operator, 

J2i=i CT W r -W; acting only on the spin and isospin, cannot change the symmetry 

properties of orbital wave functions. For Hep, therefore, the leading one-body GT 

operator cannot connect the main components of the initial and final states - 

a feature that leads to a drastic suppression of the Hep amplitude. This implies 

that the exchange-current (EXC) effects may play an exceptionally large role here. 

Werntz and Brennan (36) attempted to relate the Hep rate to the Ml matrix 

element for the Hen process, where Hen is radiative capture of a thermal neutron 

on 3 He: 3 He + n — > 4 He + 7. They assumed (a) the validity of isospin symmetry 

apart from the difference in the radial functions of the incident nucleons (proton 

for Hep and neutron for Hen); and (b) that two-body EXC terms dominated for 

both Hep and Hen and that their matrix elements could be related to each other 

via an isospin rotation. Based on the upper limit for the Hen cross section known 

at that time, Werntz and Brennan gave an upper limit for the Hep 5-factor, 3.7 x 

10~ 20 keV-b, which was about 200 times smaller than Salpeter's estimate. Later 
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Werntz and Brennan (37) refined their estimate in several respects, including the 

addition of the contributions from p-wave capture channels, and they arrived at 
an 5-factor of 8.1 x 10~ 20 keV-b. Tegner and Bargholtz (38) also attempted to 
relate Hep to Hen, but they pointed out the importance of the contributions due 
to the .D-state components in the 3 He and 4 He wave functions, and they used 
the EXC operators of the type derived by Chemtob and Rho (39); Tegner and 
Bargholtz's estimate was (4 - 25)xl0 -20 keV-b, where the spread corresponded 
to the range of the experimental values of the Hen cross section before 1983. 
This result was sharpened by Wolfs et al. (40), who measured the Hen cross 
section precisely and reported a value of (54 ±6) x 10 -20 uh; their estimate of the 
Hep 5-factor was (15.3 ±4.7) x 10~ 20 keV-b. Wervelman et al. (41) also made a 
precision measurement of the Hen cross section and obtained (55 ±3) x 10~ 20 uh, 
in good agreement with that obtained by Wolfs et al.. However, Wervelmann et 
al. predicted a Hep 5-factor of (57 ± 8) x 10 -20 keV-b. These estimates should 
be considered semi-quantitative, since even the estimates of the one-body terms 
differ wildly from model to model, and furthermore it is known (39) that the 
EXC for GT transitions should differ from that for Ml transitions. Subsequently, 
Carlson et al. (42) showed that there is a signicant cancellation between the one- 
body and two-body terms and that the use of realistic wave functions (as opposed 
to schematic wave functions employed in the previous calculations) is crucial for 
a reliable estimation of the Hep rate. These authors performed a variational 
Monte Carlo calculation and, with the use of EXC operators derived from pion- 
and /9-exchange diagrams and A-excitation diagrams, they obtained S = 1.3 x 
10~ 20 keV-b (42). Schiavilla et al. (43) performed a similar calculation but with 
the use of explicit A degree of freedom and obtained S = (1.4—3.2) x 10~ 20 keV-b. 
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The Hep calculations up to this point only considered the contribution of the s- 

wave capture channel, except in the work of Werntz and Brennan (37). Horowitz 

(44) presented a new estimate of the contribution of the p-wave capture channel, 

using schematic wave functions, and emphasized that it could be of substantial 

magnitude. We categorize Hep estimations that have appeared since 2000 as 

"recent" calculations and discuss them in the next section. 

As the above survey shows, the calculated value of the Hep 5-factor changed by 

two orders of magnitude from the original Salpeter value. Fortunately, however, 

an encouraging sign of convergence in the Hep 5-factor has been emerging over 

the past few years. This is attributable, first, to further significant progress 

along the line of work following References (42,43). Second, the application of 

effective field theory to Hep has greatly increased the reliability of the calculated 

S- factor. These latest developments are the subjects of the following sections. To 

present them coherently, we first survey the relevant theoretical frameworks in 

a somewhat general context, and then proceed to discuss the specific numerical 

results for Hep. 

3 RECENT CALCULATIONS - THEORETICAL FRAMEWORK 

3. 1 Standard Nuclear Physics Approach ( SNPA ) 

The phenomenological potential picture has been highly successful in describing 
a great variety of nuclear phenomena. In this picture a model Hamiltonian for an 
vl-nucleon system involves a phenomenological two-body potential v phen (and, if 
needed, an additional phenomenological three-body potential). Once this model 
Hamiltonian is specified, the nuclear wave function is obtained by solving the 
A-body Schrodinger equation. Recent progress in numerical techniques for this 
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type of calculation has reached such a level (45) that the wave functions of low- 
lying levels for light nuclei can now be obtained nearly without approximation. 
This achievement frees us from the "usual" nuclear physics complications that 
arise from truncation of nuclear Hilbert space down to certain model space (such 
as shell-model configurations, cluster-model trial functions, etc.). Because there 
is large freedom in choosing a possible form for the short-range part of v phen , 
one assumes a certain functional form and fixes the parameters appearing in it 
by demanding that the nucleon-nucleon (NN) scattering data and the deuteron 
properties be reproduced. There are by now several so-called high-precision phe- 
nomenological potentials that can reproduce all the existing two-nucleon data 
with normalized % 2 values close to unity. In normal circumstances, nuclear re- 
sponses to external electroweak probes are given, to good approximation, by 
one-body terms; these are also called the impulse approximation terms. To ob- 
tain higher accuracy, however, one must also consider exchange current (EXC) 
terms, which represent the contributions of nuclear responses involving two or 
more nucleons. In particular, if for some reason the impulse approximation con- 
tributions are suppressed, it becomes essential to take account of the EXC con- 
tributions. These EXC's are usually derived from one-boson exchange diagrams, 
which impose the low-energy theorems and current algebra properties on the ver- 
tices featured in the diagrams (39,46,47). A formalism based on this picture 
is referred to as the standard nuclear physics approach (SNPA), also called the 
potential model in the literature. SNPA has been used extensively to describe 
nuclear electroweak processes in light nuclei, and the generally good agreement 
between theory and experiment (45) gives a strong indication that SNPA captures 
much of the physics involved. The calculations quoted earlier (42,43) represent 
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3.2 Effective Field Theory (EFT) 

Although SNPA has been scoring undeniable successes in correlating and ex- 
plaining a vast variety of data, it is still important from a fundamental point 
of view to raise the following issues. First, since hadronic systems are governed 
by quantum chromodynamics (QCD), one should ultimately be able to relate 
nuclear phenomena to QCD, but this relation is not visible in SNPA. In par- 
ticular, whereas chiral symmetry is known to be a fundamental symmetry of 
QCD, SNPA is largely disjoint from this symmetry. Second, in SNPA, even for 
describing low-energy phenomena, we start with a "realistic" phenomenological 
potential that is tailored to encode short-range (high-momentum) and long-range 
(low-momentum) physics simultaneously. This mixing of the two different scales 
seems theoretically unsatisfactory. Third, in writing down a phenomenological 
Lagrangian for describing the nuclear interaction and nuclear responses to the 
electroweak currents, we find no clear guiding principle in SNPA — no obviously 
identifiable expansion parameter that helps us to control the possible forms of 
terms in the Lagrangian and that provides a general measure of errors in our 
calculation. To address these and related issues, a new approach based on EFT 
was proposed (48), and it has been studied with great intensity (for reviews, see 
References (49)- (50)). 

The intuitive picture of EFT is rather simple. In describing phenomena char- 
acterized by a typical energy- momentum scale Q, we expect that we need not 
include in our Lagrangian those degrees of freedom that pertain to energy- 
momentum scales much higher than Q. This expectation motivates us to intro- 
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duce a cut-off scale A that is sufficiently large compared with Q and to classify 
our fields (to be generically represented by $) into two groups: high-frequency 
fields (&high) and low- frequency fields {&iow)- By eliminating (or "integrating 
out") ^highi we arrive at an effective Lagrangian that only involves <&i ow as ex- 
plicit dynamical variables. Using the notion of path integrals, we find that the 
effective Lagrangian £ e g is related to the original Lagrangian C as 



One can show that £ e ff defined by eq.(l) inherits the symmetries (and the pat- 
terns of symmetry breaking, if there are any) of C It also follows that £ e s 
should be the sum of all possible monomials of <I>; ow and their derivatives that 
are consistent with the symmetry requirements dictated by C. Because a term 
that involves n derivatives scales like (Q/A) n , the terms in £ e ff can be organized 
into a perturbative series in which Q/A serves as an expansion parameter. The 
coefficients of terms in this expansion scheme are called the low-energy constants 
(LECs). Insofar as all the LEC's up to a specified order n can be fixed either 
from theory or from fitting to the experimental values of the relevant observables, 
£ e ff serves as a complete (and hence model-independent) Lagrangian to the given 
order of expansion. 

When EFT is applied to nuclear physics, the underlying Lagrangian is the QCD 
Lagrangian Cqcd, whereas, for the typical nuclear physics energy- momentum 
scale Q <C A x ~ 1 GeV, the effective degrees of freedom that would feature in 
the effective Lagrangian C c g are hadrons rather than the quarks and gluons. It 
is not obvious how to apply the formal definition, eq.(l), to establsh a relation 
between Cqcd and £ e ff written in terms of the hadrons, since the hadrons cannot 
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be simply identified with the low-frequency field, <&i ow , in the original Lagrangian. 

At present, the best one can do is to resort to symmetry considerations and the 

above-mentioned expansion scheme. Here chiral symmetry plays an important 

role. It is known that chiral symmetry is spontaneously broken, generating the 

pions as Nambu-Goldstone bosons. This feature can be incorporated by assigning 

suitable chiral transformation properties to the Goldstone bosons and writing 

down all possible chiral-invariant terms up to a specified chiral order (see e.g. 

Reference (51)). The above consideration presupposes exact chiral symmetry 

in £qcd- In reality, £qcd contains small but finite quark mass terms, which 

explicitly violate chiral symmetry and lead to a non-vanishing value of the pion 

mass m w . Again, there is a well-defined method to determine what terms are 

needed in the Goldstone boson sector to represent the effect of explicit chiral 

symmetry breaking (51). We can then establish a counting rule, called chiral 

counting, which allows us to classify the relative importance of a term in C e g 

and a given Feynman diagram according to the number of powers in Q/A and 

m w /A. These considerations lead to an EFT called chiral perturbation theory 

( X PT) (52,53). 

The successes of xPT in the meson sector are well known (see, e.g., Reference 
(49)). A problem we encounter in extending xPT to the nucleon sector is that, 
because the nucleon mass is comparable to the cut-off scale A x , a simple ap- 
plication of expansion in Q/A does not work. This problem can be circumvented 
by employing heavy-baryon chiral perturbation theory (HB^PT) (54). HB^PT 
has been applied with great success to the one- nucleon sector (see e.g. Refer- 
ence (49)). HB^PT, however, cannot be applied in a straightforward manner 
to nuclear systems, because nuclei involve very low-lying excited states, and the 
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existence of this small energy scale spoils the original counting rule (48). 



3.3 Nuclear EFT in the Weinberg Scheme - K-Counting 

Weinberg proposed to avoid this difficulty by classifying Feynman diagrams into 
two groups, irreducible and reducible diagrams (48). Irreducible diagrams are 
those in which every intermediate state has at least one meson in flight; all 
others are classified as reducible diagrams. We then apply the above-mentioned 
chiral counting rules only to irreducible diagrams. The contribution of all the 
two-body irreducible diagrams (up to a specified chiral order) is treated as an 
effective potential acting on nuclear wave functions. Meanwhile, the contributions 
of reducible diagrams can be incorporated by solving the Schrodinger equation. 
This two-step procedure may be referred to as nuclear xPT or, to be more specific, 
nuclear xPT in the Weinberg scheme. 

To apply nuclear %PT to a process that involves (an) external current (s), we 
derive a nuclear transition operator T by evaluating the complete set of irre- 
ducible diagrams (up to a given chiral order v) that involve the relevant external 
current (s). For consistency in chiral counting, the nuclear matrix element of T 
must be calculated with the use of nuclear wave functions that are governed by 
nuclear interactions that represent all the irreducible A-nucleon diagrams up to 
z/-th order. If this program could be carried out exactly, it would constitute an 
ab initio calculation. The unambiguous classification of transition operators ac- 
cording to their chiral orders is a great advantage of EFT, which is missing in 
SNPA. 
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3.4 Nuclear EFT in the KSW Scheme - Q-Counting 

An alternative form of nuclear EFT is based on the power divergence subtraction 
(PDS) scheme. The PDS scheme proposed by Kaplan, Savage and Wise (KSW) 
in their seminal papers (55) uses a counting scheme (often called Q-counting) 
that differs from the Weinberg scheme. Although a great number of important 
investigations have used the PDS scheme (for a review, see e.g. Reference (56)), 
here we are primarily concerned with the Weinberg scheme. The reason is that 
the PDS scheme has so far been used chiefly only for the two-nucleon systems 
(see, however, References (57-59)), and so at present it is less directly connected 
with the Hep process than the Weinberg scheme is; no Hep calculations based 
on the PDS scheme exist in the literature at the time of this writing. 

3.5 Hybrid EFT 

The preceding subsections emphasize the formal merits of nuclear EFT. In actual 
calculations, however, it is still a major challenge to generate, strictly within the 
EFT framework, nuclear wave functions whose accuracy is comparable to that 
of SNPA wave functions (see section 7, however). A pragmatic solution to this 
problem is to use wave functions obtained in SNPA; we refer to this eclectic 
approach as hybrid EFT (60-62). Since the AW interactions that generate SNPA 
wave functions accurately reproduce the two-nucleon data, the use of hybrid-EFT 
is almost equivalent to using the empirical data themselves to control the initial 
and final nuclear wave functions, insofar as the off-shell problem (see below) 
and the contributions of three-body (and higher-body) interactions are properly 
addressed. 
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3.6 EFT* or MEEFT 

Hybrid EFT can be used for light complex nuclei (A = 3, 4, ...) with essentially 
the same accuracy and ease as for the A=2 system. We should emphasize in this 
connection that, in A-nucleon systems (A > 3), the contributions of transition 
operators involving three or more nucleons are intrinsically suppressed according 
to chiral counting, and hence, up to a certain chiral order, a transition operator 
in an ^4-nucleon system consists of the same EFT-based one-body and two-body 
terms that are used for the two-nucleon system. 

As mentioned above, the chiral Lagrangian is definite only when the values of all 
the relevant LECs are fixed, but there may be cases where this condition cannot be 
readily met. Suppose that a two-body EXC operator under consideration contains 
an LEC (call it k) that cannot be determined with the use of A=2 data alone. It 
is possible that an observable (call it fi) in a A-body system [A > 3) is sensitive 
to k and that the experimental value of Q is known with sufficient accuracy. 
Then we can determine k by calculating the hybrid-EFT matrix element that 
corresponds to Q and adjusting k to reproduce the empirical value of Q. Once k 
is fixed this way, we can make predictions for any other observables for any other 
nuclear systems that are controlled by the same transition operators. Hybrid 
EFT used in this manner is referred to as EFT* or as MEEFT (more effective 
EFT). 

The effective Lagrangian C e ff is, by construction, valid only below the spec- 
ified cutoff scale A. This basic constraint should be respected in nuclear EFT 
calculations as well. One way to implement this constraint is to introduce a 
momentum- cut off A for the two-nucleon relative momentum. The sensitivity of 
the results on the choice of A is expected to serve as a measure of uncertainties 
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in the calculational framework. 
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4 NUMERICAL RESULTS 

4-1 Hep Calculation Based on SNPA 

Marcucci et al. have recently carried out a highly elaborate SNPA calculation of 
the Hep rate (34,63). The treatment of the four-body wave functions was im- 
proved with the use of the correlated- hyper spherical-harmonics method (64). The 
strength of the dominant EXC contribution due to a A-excitation diagram was 
adjusted to reproduce the experimental value of r^ ritium , the tritium beta-decay 
rate: rjf tium (exp) = 12.32 ± 0.03 yr (65). This type of empirical normaliza- 
tion, first introduced in Reference (42), is expected to reduce significantly the 
model dependence of the calculated Hep rate. Furthermore, the contribution 
of the initial p-wave channel was included. The resulting Hep 5-factor at 10 
keV (close to the Gamow peak) is S = (10.1 ± 0.6) x 10~ 20 keV-b; this is the 
value used by Bahcall et al. (BP00) (3). The corresponding threshold value is 

5 = 9.64 x 10~ 20 keV-b. 

4-2 Hep Calculation Based on EFT* 

Park et al. (66, 67) carried out an EFT* calculation of the Hep rate up to next- 
to-next-to-next-to-leading order in chiral counting. To this order, there appears 
in two-body terms one LEC that at present cannot be determined from data 
belonging to the A=2 systems. This unknown LEC, denoted by d,R in Reference 
(62), parametrizes the strength of a contact-type four-nucleon axial-current cou- 
pling. Park et al. noted that dn also appears as the only unknown parameter 
in the calculation of ri ntmm . They determined dn from the experimental value 
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of r^ ntmm (65). With the value of dn determined this way, Park et al. made a 

parameter- free EFT* calculation of the Hep rate (66,67). The result for the 
threshold S-factor is S = (8.6 ± 1.3) x 10~ 20 keV-b, where the error spans the 
range of the A dependence. This EFT* result supports the SNPA results in 
References (34,63). It is reasonable to expect that, if the effects of neglected 
higher-order terms and deviations from the framework of EFT* itself are sizable, 
they would cause the significant A dependence in the calculated S. This consid- 
eration led Park et al. (66,67) to adopt the A dependence (~15 % variation) of 
Sjiep as a measure of uncertainties in their EFT* calculation of Suep- The above- 
mentioned large cancellation between the one-body and two-body contributions 
in Hep amplifies the A dependence of Suep up to ~15 %, but this dependence is 
still reasonably small. (Below we discuss the pp-fusion reaction, where there is 
no such cancellation and the A dependence is found to be extremely small.) In 
Fig.2, the 15 % errors assigned to 4> v {hep) reflects the uncertainty in the EFT* 
estimation of Snep in References (66,67); this is the first time that 4> u (hep) is 
presented with an error estimate attached. 

4.3 Off- Shell Problem 

The use of hybrid EFT may bring in a certain degree of model dependence due to 
off-shell effects, because the phenomenological NN interactions are constrained 
only by the on-shell two-nucleon observables. This off-shell effect, however, is 
expected to be small for the reactions under consideration, since they involve low 
momentum transfers and hence are not extremely sensitive to the short-range 
behavior of the nuclear wave functions. One way to quantify this expectation is 
to compare a two-nucleon relative wave function generated by the phenomeno- 
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logical potential with that generated by an EFT-motivated potential. Phillips 

and Cohen (68) made such a comparison in their analysis of the one-body oper- 
ators responsible for electron-deuteron Compton scattering. They showed that 
a hybrid EFT should work well up to momentum transfer 700 MeV. A similar 
conclusion is expected to hold for a two-body operator, so long as its radial be- 
havior is duly "smeared-out" reflecting a finite momentum cutoff. Thus, EFT* 
as applied to low energy phenomena is expected to be practically free from the 
off-shell ambiguities. 

Another indication of the stability of the EFT* results comes from the recently 
developed idea of the "low-momentum nuclear potential". Let us recall that 
a "realistic phenomenological" potential v phen is determined by fitting to the 
two-nucleon data up to the pion production threshold energy. So, physically, 
v phen ghQuid re side in a momentum regime below a certain cutoff A c . In the 
conventional treatment, however, the existence of this cutoff scale is ignored. 
Bogner et al. (69) proposed to construct an "effective low-momentum" potential, 
Vlow-k, by integrating out from v phen the momentum components higher than 
A c . They calculated V/ ow _fcS corresponding to a number of well-established NN 
potentials. Remarkably, the resulting Vj ow _fcS were found to lead to identical 
half-off-shell T-matrices for all the cases studied, even though the ways short- 
range physics is encoded in these v phen s are quite diverse. This implies that the 
Vi ow -kS are free from the off-shell ambiguities, and therefore the use of Vi ow -kS is 
essentially equivalent to employing an EFT-based AW potential. The fact that 
EFT* calculations by design contain a momentum-cutoff regulator essentially 
ensures that an electroweak transition matrix element calculated in EFT* is only 
sensitive to those half-off-shell T-matrices that are controlled by Vi ow -k, and 
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therefore the EFT* results reported by Park et al. (67) are expected to be 

essentially free from the off-shell ambiguities. 

Furthermore, because correlating the observables in neighboring nuclei (as was 
done here) is likely to serve as an additional renormalization, the possible ef- 
fects of higher-chiral-order terms and/or off-shell ambiguities can be significantly 
suppressed by the use of EFT* 

5 COMPARISON WITH EXPERIMENTAL DATA 

In Super-Kamiokande experiments, information on the solar neutrino spectrum 
4> u was obtained by detecting the recoil electron in the reaction, v + e~ — > v + e~, 
for -E r ccoii > 6.5 MeV (29), and for -E reco ii > 5 MeV (14). As mentioned, 4> u in this 
range is governed by the 8 B neutrinos mixed with a tiny number of hep neutrinos. 
Reference (29) presents the data in 15 bins between 6.5 MeV and 14 MeV and 
one higher-energy bin covering 14 - 20 MeV. The three highest bins showed a 
larger number of events than expected from the then most popular neutrino 
oscillation parameters. Bahcall and Krastev (33) analyzed these data in detail; 
they considered various neutrino oscillation scenarios and treated the Hep S- 
factor, Sjjep, as an adjustable parameter. The philosophy behind this treatment 
was that, although the "1998 standard value" of Suep adopted in (70) came from 
an elaborate SNPA calculation (43), a first-principle physics argument was still 
needed to exclude the possibility that Suep might exceed, e.g., 10 times this value. 
Introducing the enhancement factor a defined by a = S# ep /(2.3 x 10~ 20 keV-b), 
where the denominator is the central value of the 1998 standard Sn e p, Bahcall and 
Krastev reported that, by allowing a to be larger than 20, one could significantly 
improve global fits to all the then available solar neutrino data for every neutrino 
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oscillation scenario studied. This result triggered renewed theoretical efforts by 

Marcucci et al. (34, 63) and Park et al. (66, 67) to determine Suep with higher 
accuracy. An improved SNPA calculation by Marcucci et al. (34,63) gave a value 
of Suep 4.4 times larger than the "1998 standard value"; the value of (f) u (hep) SSM 
that appears in BPOO (3) is based on Marcucci et al.'s Sn e p- It is to be noted that 
the authors of (3) avoided giving an estimate of total uncertainty in <p u {hep) , 
citing the unique subtlety involved in the calculation of S 'Hep- 

According to the more recent Super-Kamiokande results (14) covering -E reC oii > 
5 MeV, the observed shape of <p v is consistent with an undistorted 8 B neutrino 
spectrum shape; a \ 2 fit of the overall spectrum shape resulted in x 2 /d.o.f. = 
19.1/18, without any Hep neutrino admixture. This is not very surprising since 
4> v {hep) is typically three orders of magnitude smaller than <f> u ( 8 B). Meanwhile, 
by assuming that 1.3 ± 2.0 events registered in the recoil electron energy bin, 
-E'recoil = 18 — 21 MeV, are due to the Hep neutrinos, the 90 % confidence level 
upper limit of 4> v (hep) was determined to be 40 x 10 3 cm _2 s _1 (14). This upper 
limit is 4.3 times the BPOO prediction for the no neutrino-oscillation assumption; 
thus the experimental (f) u (hep) is consistent with the theoretical value. However, 
BPOO stressed that the significance of this agreement is limited because one can- 
not assign an estimate of uncertainty to the theoretical value of Snep used in 
BPOO. 

This drawback can be greatly mitigated with the use of the EFT* calculation 
of Suep by Park et al. (66,67), which gives Sh&p with a well-controlled error 
estimate (in the sense explained earlier). The use of the central value of Sjiep 
obtained by Park et al. (67) would slightly lower 4> u (hep) but would keep its 
upper end compatible with <f> u (hep) SSM in BPOO (3). Thus, the statement that 
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the upper limit of the experimental (j) u {hep) is ~4 times the central value of the 

SSM prediction remains valid. However, with the use of Suep obtained in the 

EFT* calculation, the SSM prediction is controlled within ~15 % precision, and 

this fact is expected to be valuable in future analyses of experiments concerning 

Hep neutrinos. Coraddu et al. (71) pointed out that, with the precision of Shep 

achieved in Reference (67), we may be able to use the 4> u (hep) data to study the 

possible nonstandard behavior of the solar core plasma. 

6 RELATED TOPICS 

We have so far concentrated on the calculations of Hep. To further clarify the 
key aspects involved in these calculations, we now discuss the related problems 
of evaluating the cross sections for the pp fusion reaction (p + p — > d + e + + v e ), 
the neutrino-deuteron reactions, and the Hen reaction. 

6. 1 The pp Fusion Process 

The latest SNPA estimation of the pp fusion rate was carried out by Schiavilla et 
al. (72), using the EXC operator whose strength was normalized to fit r^ ritlum . 
Park et al. (67, 73) performed an EFT* calculation of the pp fusion rate, us- 
ing exactly the same method employed for the Hep calculation. The result is 
S pp = 3.94 x (1 ± 0.005) x 10~ 25 MeVb. This EFT* result supports the value of 
S pp obtained in SNPA (72). It has been found that S pp in the EFT* calculation 
varies by only ~0.1 % against changes in A, and this feature can be regarded 
as typical for unsuppressed transitions such as pp fusion. Thus, to the extent 
that the A dependence serves as a reasonable measure of theoretical uncertain- 
ties (as discussed above), the EFT* prediction of Spp can be considered highly 
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robust. The 0.5 % error in the above S pp is dominated by the uncertainty in the 

experimental value of r^ ritmm . 

The PDS scheme also was used to estimate the pp fusion rate (74). Taking 
advantage of very low energy and momentum involved in this reaction, Kong 
and Ravendal used "nucleon-only" EFT (without pions). To the order they con- 
sidered, there appears only a single LEC, denoted by Lia, and there have been 
attempts to constrain its value using observables in the two-nucleon systems (75). 

6.2 Neutrino- deuteron reactions 

The u-d cross sections for E v < 20 MeV are very important in connection with 
the SNO experiments. Nakamura et al. performed a detailed SNPA calculation 
of the v-d cross sections o(yd) (76,77), and Butler et al. (78) carried out an EFT 
calculation of a(vd), using the PDS scheme (55). The EFT results (78) agree 
with the SNPA results (77), if the above-mentioned unknown LEC, Lia, involved 
in the PDS scheme is suitably adjusted. The optimal value, Lia = 5.6 fm 3 , found 
by Butler et al. (78) is consistent with the order of magnitude of Lia expected 
from the naturalness argument (based on a dimensional analysis), |Lia| < 6fm 3 . 
Even though it is reassuring that a(ud) calculated in SNPA and EFT agree 
with each other (in the above-explained sense), it is desirable to carry out an 
EFT calculation that is free from any adjustable LEC. EFT* allows us to carry 
out an EFT-controlled parameter-free calculation of a(vd), and Ando et al. (79) 
performed such a calculation. The a(vd)s they obtained (79) are found to agree 
within 1% with a(udys obtained in SNPA (77). Although it is in principle possible 
to calculate the tritium beta decay rate in the PDS scheme and use ^^^(exp) 
to determine Lia, this program has yet to be carried out. 
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6.3 Hen process 

In order to gauge the validity of a calculational method used for Hep, it is ex- 
tremely useful to study as a test case the Hen process, 3 He + n — > 4 He + 7. This 
is because Hep and Hen involve similar kinematics and share the characteristic 
that the contribution of the normally dominant one-body transition operator is 
highly suppressed owing to the symmetry properties of the wave functions. An 
EFT* calculation of Hen has been carried out by Song and Park (80), and the 
calculated cross section, a = (60.1 ±3.2 ±1.0) fib, is in reasonable agreement with 
the experimental values, (54 ± 6) ub (40), and (55 ± 3) fib (41). This agreement 
supports the validity of the EFT* approach in general and the EFT* calculation 
of Hep in particular, even though there is room for improvements in the j>- 3 He 
continuum wave function used in Reference (80) . For the earlier Hen calculations 
based on SNPA, see References (43,81). It is hoped that there will be further 
investigations of Hen in both SNPA and EFT*. 

7 SUMMARY AND OUTLOOKS 

As exemplified above, low-energy electroweak processes in light nuclei play im- 
portant roles in astrophysics, and a recently developed EFT-based formalism, 
called EFT* or MEEFT (more effective EFT), can be used profitably to calcu- 
late the cross sections of these processes with high precision. We have discussed 
here Hep and a few closely related reactions, but this new method is expected to 
prove useful for other low-energy electroweak processes as well. The numerical 
results obtained in EFT* generally support those obtained in the conventional 
SNPA, if the strength of the two-body current is controlled by the empirical value 
of an appropriate observable. It is to be stressed that EFT* allows us to make 
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systematic error estimation of the calculated cross sections, a feature that is not 

readily obtainable in SNPA. 

From a formal point of view, one could hope to improve EFT* by employing 
nuclear wave functions determined in an EFT formalism itself instead of phe- 
nomenological wave functions obtained in SNPA. In regard to observables that 
do not involve external currents, there has been great progress in building a for- 
mally consistent EFT approach applicable to complex nuclei (58,82). It will be 
highly informative to apply this type of formalism to electroweak processes and 
compare the results with those of EFT*. 
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p + p^ 2 H + e+ + u e (99.75%) p + e~ + p -> 2 H + u e (0.25%) 



I 

2 H + p^ 3 He + 7 

3 He + 3 He -> a + 2p (86%) 3 He + a -> 7 Be + 7 (14%) 3 He + p -> a + e+ + ^ e 



7 Be + e" -> 7 Li + ^ e (99.85%) 7 Be + p -> 8 B + 7 (0.15%) 

I I 
7 Li + p -► 2a 8 B -> 8 Be* + e+ + u e 

I 

8 Be* -> 2a 



Figure 1: Solar thermonuclear reactions in the pp-chain and their branching 
ratios. The Hep branching ratio is of the order of 0.01% or less. 
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Figure 2: Solar neutrino spectrum 4> v versus the neutrino energy E u . The neu- 
trino fluxes from continuum sources are given in units of cm -2 s -1 MeV -1 , and 
the line fluxes in units of cm~ 2 s _1 . The arrows at the top indicate the ranges of 
E v covered by the experiments mentioned in the text. 



